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Fractional diffusion equations are widely used to describe anomalous 
diffusion processes where the characteristic displacement scales as a power 
of time. For processes lacking such scaling the corresponding description 
may be given by distributed-order equations. In the present paper we con- 
sider different forms of distributed-order fractional kinetic equations and 
investigate the effects described by different classes of such equations. In 
particular, the equations describing accelerating and decelerating subdiffu- 
sion, as well as the those describing accelerating and decelerating superdif- 
fusion are presented. 

PACS 05.40.-hj; 02.50.-r 

1. Introduction 

The diffusion equation proposed by Adolf Fick almost 150 years ago 
is a partial differential equation of parabolic type, with the first temporal 
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derivative on its l.h.s. and the second spatial derivative on its r.h.s.. The 
corresponding equation can be put down both for the particles' concentra- 
tion in a diffusing cloud and for a probability of a single particle's position: 

^^p{x,t) =K-^p{x,t) (1.1) 

The mean squared particle's displacement from its initial position given 
by the solution of this equation grows linearly in time: (x^(t)) = 2Kt. This 
scaling behavior follows immediately from the structure of Pick's equation, 
being second order in spatial coordinate and first order in time: Changing 
the spatial scale by a factor of 2 corresponds to changing the time scale by 
a factor of 4. 

In complex systems this kind of behavior is often violated, being replaced 
by an anomalous diffusion relationship, 

{xHt))^t^, (1.2) 

with P ^ 1, or (in absence of the second moment) by related forms where the 
lower, not necessarily integer, moment of the distribution scales as a func- 
tion of time. The continuous description of anomalous diffusion is given 
by generalizations of Pick's scheme based on fractional derivatives PP. The 
fractional generalizations of the diffusion equation may have either a corre- 
sponding fractional derivative instead of the whole-number one, or an ad- 
ditional fractional derivative on the "wrong" side of the equation. In what 
follows we refer to these two possibilities as to "normal" and "modified" 
forms of a fractional diffusion equation. An example of the "normal" form 
is an equation for superdiffusion with the Riesz fractional spatial deriva- 
tive instead of the second derivative on the r.h.s. 2 . The situation with 
an additional derivative on the "wrong" side is exemplified by the stan- 
dard fractional diffusion equation for subdiffusion [Sj. We note that the two 
forms, i.e. the form with a Caputo derivative on the l.h.s. for subdiffusive 
processes, and a form with an additional spatial derivative on the l.h.s. for 
a super diffusive process, are equivalent to the commonly used ones. This 
equivalence will be discussed in detail in the Section 2 of the present work, 
introducing the corresponding fractional operators. 

Many physical processes, however, lack power-law scaling, Ea. (|1.2|l . over 
the whole time-domain. These processes can not be characterized by a single 
scaling exponent p. Examples of processes lacking scaling include several 
cases of decelerating subdiffusion (e.g. the Sinai diffusion j^) and decelerat- 
ing superdiffusion (as exemplified by truncated Levy- flights [3]). Such pro- 
cesses can be described by derivatives of distributed order 017], introduced 
by Caputo |S|. The equivalence between the different forms of fractional dif- 
fusion equations is lost in this case: Different forms of the distributed-order 
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fractional equations describe different situations. Thus, the equations with 
the distributed-order derivatives on the "proper" side describe processes 
getting more anomalous in course of the time (accelerating superdiffusion 
and decelerating subdiffusion) while the equations with the additional 
distributed-order on the "wrong" side describe the situations getting less 
anomalous (decelerating superdiffusion and accelerating subdiffusion). As 
an example a special model with two fractional derivatives of different orders 
is used throughout the work. 

In Section 2 and 3 we discuss the four forms of fractional diffusion equa- 
tions for the processes showing scaling behavior. These forms correspond 
to possible permutations of fractional temporal/spatial derivative on the 
"proper" /"wrong" side of the equation. Three of the four forms discussed 
below are well-known. The fourth one, to our best knowledge, was not 
previously discussed. We first turn to the forms pertinent to temporal frac- 
tional equations, i.e. to subdiffusion, and then to the superdiffusive case. 
Sections 4 and 5 are devoted to the distributed-order generalizations of the 
corresponding equations. The results are summarized in Section 6. 



2. Two forms of time fractional diffusion equations 

2.1. Riemann-Liouville form 

The time fractional diffusion equation (TFDE) in the Riemann-Liouville 
form (RL - form), which in our terminology corresponds to a "modified" 
form of the fractional diffusion equation, reads 

-pix, t) = oD]-^^p{x, t), (2.1) 

p{x,t = 0) = Six), < /3 < 1. 

Here Kfj is a positive constant, [K/^] = cm^/sec^, and qD^ is the Riemann- 
Liouville fractional derivative on the right semi-axis, which, for a "suffi- 
ciently well-behaved" function (j){t) is defined as follows: 



d ^ ^ rt 

di" ^~ T{l-fi)dt Jq"" (t-rY' 



oD^cP = -J'-^4> = / ^^77^' < < 1, (2.2) 



where J"(j){t) = JodT{t - t)°'-'^ <j){t) , t > 0, a E M+ is the Riemann- 
Liouville fractional integral of the order a. In what follows, we omit the 
subscript "0", for brevity. 

Applying the Laplace transform, 

POO 

4>is)^L{cl>{t)}= / dte-^V(t), 
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and Fourier-transform, 



/oo 
-oo 



in succession, and using the Laplace transform of the Riemann-Liouvihe 
derivative ()2.2() . 

L{D['<p{t)] = s^ks) (2.3) 

we get from Ea. H2.1() the form of the Laplace-transformed characteristic 
function f{k,t) of the distribution p{x,t): 

2.2. Caputo form 

The TFDE in the Caputo form (C - form, corresponding to a "normal" 
form of the fractional diffusion equation) is written as follows: 

Q^p{x, t) = Kfs-^p{x, t), (2.5) 

p{x,0) = 5{x), where Kp is the same constant as in Ea. (|2.1() . and the time 
fractional derivative of order /3, < /? < 1 is understood in the Caputo 
sense [TT] . 



Here, the sequence of temporal integration and differentiation is reversed 
with respect to a Riemann-Liouville operator. Recalling the Laplace trans- 
form of the Caputo derivative, 

L[^]=s^ks)-s^-'m, (2.7) 

< < 1, and making the Fourier-Laplace transform of Ea. (|2.5|) we again 
arrive at Ea. H2.4|) . Thus, both forms of TFDE are equivalent. 

In the literature, the third form of TFDE is mentioned, this one also 
using the Riemann - Liouville derivative jl2j : 

D^^p{x, t) = Kp^p{x, t) + -l^p(x, 0), (2.8) 
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< /3 < 1. The equivalence of the RL- and this third forms can be easily 
shown by applying dt' ... to both sides of Ea. H2.1() and using the Riemann 

- Liouville derivative of a constant, -Of 1 = i~^/r(l — The equivalence of 
the "normal" (C-) and the third forms can also be shown easily, if one uses 
the relation between the Riemann-Liouville and Caputo derivatives which 
can be obtained straightforwardly: 

with < /3 < 1. Then, starting, e.g., from the C - form of TFDE, Eq. (|2.5|) . 
and using Ea. (|2.9jl we immediately arrive at Ea. (|2.8|) . 

In this paper we are interested in the mean squared displacement (MSD) 
given by 

(-^W> = [ \ ■ (2.10) 



k=0. 



From Ea. H2.4|l we get: 



(^^t))=L-^{2K,.-'-^}^^^. (2.11) 

3. Two forms of space fractional diffusion equations 
3.1. The "normal" form 
The "normal" form of space fractional diffusion equation reads as: 

dp _ . . 

p{x,t = 0) = 5{x), < a < 2, where Ka is a positive constant, [Ka] = 
cm" /sec, and the Riesz fractional derivative d°^/d\x\°' (we adopt here the 
notation introduced in is defined for a "sufficiently well-behaved" func- 
tion /(x) through the Liouville - Weil derivatives jlUj : 



^/W- (3.2) 

^1^1 1 -iiHf{x) a = l 



where are the left - and right side Liouville - Weil derivatives, 

1 ^2 r (l){£)di 



D 



a 



r(2 - a) dx^ J^^ (^ - 



-'^ - w^)d^i {X - ir-^ ^^-^^ 
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for 0<a<2, a^l (for a = 1 D\. = ±d/dx), and H is the Hilbert 
transform operator, 



In Fourier space the operators of fractional derivatives have a simple form: 

/oo 
dxexp{ikx)D1(P= {^ik)'^<l){k), (3.4) 
-oo 



where 



Since 



am 



(ipiA:)" = |/c|"exp ( =F— sign/c ) . 



<^ {H (p) = isign k(p{k) (3.5) 



then, with the use of Eas. (|3.2|) - (|3.5jl we get the expression, which is valid 
for the Fourier transform of the Riesz fractional derivative for all values of 
a: 

^(^-^^=-\krm- (3.6) 

Applying the Fourier-transform to Eq. (|3.1|) , and noting Eq. (j3.6j) , we get the 
characteristic function for the PDF of Levy flights, 

fik,t)=expi-D\krt). (3.7) 

3.2. The "modified" form 

Let us turn to the fractional equation for superdiffusion with the addi- 
tional spatial derivative on its l.h.s. 

Q2-a Qp Q2 

-Ko^T^P^ (3-8) 



where Ka is the same as in Ea. ()3.1|) . Note the minus sign in Ea. (|3.8() . This 
sign gets clear when turning to a Fourier representation: applying Fourier 
transform and using Ea. ()3.6() . we arrive at the characteristic function 1)3. 7|) . 
Thus, both forms of the space fractional diffusion equation, Eas. H3.1() and 
(|3.8j) . are equivalent. 

Since the mean square displacement diverges for Levy flights, their 
anomalous nature can be characterized by a typical displacement 5x of the 
diffusing particle, 

5x oc q<a, (3.9) 
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which differs, of course, from the MSD discussed in Eq. (|2.1()j) . To get the 
q-th. moment we use the following expression |13j : 



oo 



dxl") = -r(l + g)sin^ / dk{l-Ref{k,t))k-'^-\ (3.10) 



vr ' 2 

Inserting Ea. (|3.7|) into Ea. ()3.1U() and changing to a new variable, = K^k'^t, 
we obtain 

Jo Oi Jo q \ aJ 

(3.11) 

dxl") = C7(g,a)(E:«t)«/", g<a. (3.12) 



and, thus 



Here 

C{q, a) = — sin (^) r(l + g)r (l - ^) • (3.13) 
nq \ 2 / \ a) 

Note that for a = g = 2 Eqs.(jSIIll), ^^^>^ give 

(x^) = 2iv:t. 

4. Distributed-order time fractional diffusion equations 

J^.l. Distributed- order time fractional diffusion equation in the RL form 

The fractional diffusion equation with a distributed-order Riemann- 
Liouville derivative reads: 



p{x,0) = 6{x), where K{f5) = Kt^ ^, [K] = cm^/sec, [r] = sec, w{P) is 
a dimensionless non- negative function, which should fulfill Jq dj3w{(3) = 1. 
If we set w{f3) = 5{P — Pq), < (3o < 1, then we arrive at time fractional 
diffusion equation in the RL form, see Ea. (|2.1|) . where Kp = Kt'^~^'-\ 

We now prove that the solution of Ea. (|4.1|) is a PDF. The derivation 
here follows the method used in Its aim is to show that the random 
process whose PDF obeys Ea. (|4.1|) is subordinated to the Wiener process. 
Making a Fourier-Laplace transform of Eq. 1)4. 1|) and using Eq. (|2.3() we get 

ms)= (j-y.,^ y (4-2) 
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where ^ 

Irl{st)= [ d(3{sT)-Pw{(3). (4.3) 

JO 

We rewrite Ea. H4.2|) as follows: 

-| /"OO poo 

f(k, s) = — — / dnexp [-U [I^l + k^Kr)] = / due-''^''^^GRL{u, s) 
sIrl Jo Jo 

(4.4) 

where 

1 r M 1 

GRLiu,s) = — — ^exp -- — — (4.5) 

is a Laplace transform of a function whose properties will be specified below. 
Now, with the help of Eas. ()4.4() and ()4.5() the PDF p{x, t) can be written as 

Pix,t) = r ^e-^'^ [ ^e^' rdue--''^^GRUu,s) = 
J-oo 2vr Jsr 27rz Jq 

oo ^—x^/4:uKt 

du GRL{u,t) (4.6) 

\'AttuKt 

where Br denotes the Bromwich integration contour. 

In order to prove the positivity of p{x, t) we demonstrate that the func- 
tion GRL{u,t) is the PDF providing the subordination transformation from 
time scale t to time scale u, that is, GRL{u,t) is positive and normalized 
with respect to u for any t. At first we demonstrate normalization, 



oo 




duGRLiu, t) =L-U / du— }=L-U-} = 1. (4.7) 





Now, to prove the positivity of function GRL{u,t), according to the Bern- 
stein theorem ^^I) it is enough to show that Grl{u,s), Ea. ()4.5|) . is com- 
pletely monotonic as a function of s on positive real axis, i.e., it is positive 
and the signs of its derivatives alternate. We do this for the special case 
when 

w{P) = Bi6ip - pi) + B2d{P - 132) (4.8) 

with < /?! < /32 < 1, 5i > 0, ^2 > 0, Bi + B2 = 1. This choice allows 
us to show in a simple way the property of the diffusive behavior governed 
by distributed-order diffusion equations. This is why this case is repeatedly 
discussed in our article. 

Inserting Ea. H4.8|) into Ea. (|4.3() we have 

Irl{st) = bis-^' + b2S-^^ (4.9) 
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where bi = -Bi/r^i, 62 = B2/t^^. From Eq. (|I3|) it follows that Grl{u, s) is 
a product of two functions: Grl{u, s) = (j)i{s)(j)2{u, s) with 



02 = exp 



sIrl{st)' 



u 



Irl{st) 



(4.10) 



We will prove that both functions and 02 are completely monotonic, and 
therefore, as a product, Grl{u^ s) is completely monotonic too. 

Let us start from 0i. It can be rewritten as 0i = 0(/i(s)), where 4'{y) = 
1/y is a completely monotonic function, and h{s) = bis^~^^ + b2S^~^^ is a 
positive function with a completely monotonic derivative (this is evident by 
direct inspection). Therefore, </>i(s) is completely monotonic according to 
Criterion 2 in [13, Chapter XIII, §4. 

Now turn to 02- Again, it can be rewritten as 4>2 = (p{^l^{s)), where 
(t>{y) = exp(—uy) is a completely monotonic function. According to the 
same Criterion, it is enough to show that the positive function 

1 1 T 

ilj(s) = t = t = I t f4 11) 

Ib.l{st) 615-/31+625-/32 62 1 + ^s/32-/3i ^ 

possesses completely monotonic derivative. By differentiating £0. 1)4.11(1 we 
get 

^'(s) = ^ gfe &(/92-A)5^^-^^-^ _ 

62 l + ^s/52-/3i l + ks/32-/3i ■ ^ • ^ 

Denoting ^ = /32 — < ^ < 1, we note that this function is a product of 
three completely monotonic functions, and thus is itself completely mono- 
tone, too. Indeed, (i) the first function is s^^~^hemg a negative power of s; 
(ii) the second one is completely monotonic since it has the form g{h{s)), 
with g{y) = 1/(1 + y) being completely monotone, and h{^) = is a posi- 
tive function with a completely monotone derivative; (iii) the last function 
has the same form with g(y) = (1 + cy)/{l + dy), where < c < d with 
c = 6i/3i/62/32) d = 61/62- The function g{y) is completely monotonic since 
its n-th derivative (obtained using the Leibnitz rule) has the form 

Jn)(._ i-^r~Hc-d)nld-^ 
9 W (l + dy)"+i 

and the signs of the successive derivatives alternate. Thus, we have proved 
that ip'{s), Ea. (|4.12|) . is a completely monotonic function. Therefore the 
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exponential (j)2-, Eq. (|4.1()j) . is a completely monotonic function, too, and the 
whole function, Grl{u, s), is completely monotone as a product of two com- 
pletely monotonic functions, (pi and (/)2. Therefore, the function Grl{u-, t) is 
a PDF and, according to Ea. ()4.6() . the function t) is a PDF, too. Thus, 
we have proved that the solution of the distributed order time fractional 
diffusion equation in the Riemann-Liouville form is a PDF. 
We are interested in the MSD which is given by 

For the particular case (|4.8|) we get, by inserting Eq. (|4.9|) into Eq. (|4.13|) and 
making an inverse Laplace transform, 

where ki = BiKt and K2 = B2KT. Since < /3i < /32 < 1, at small 
times the first term in the r.h.s. of Eg. 1)4. 14(1 prevails, whereas at large 
times the second one dominates. Thus, the overall behavior corresponds to 
accelerating subdiffusion. 

4-. 2. Distributed- order time fractional diffusion equation in the C - form 

The distributed-order time fractional diffusion equation in the "normal" 
form can be written as 

/'*.^-..(«§^A-g. (4.15) 

p{x, 0) = S{x), where r is a positive constant representing some characteris- 
tic time of the problem (vide infra), [r] = sec, K is the diffusion coefficient, 
[K] = cm^/sec, w{f3) is a dimensionless non-negative function, and the time 
fractional derivative of order /3 , < /3 < 1 is understood in the Caputo 
sense. Eg. 1)2. 6(1 . Note the difference between Eg. 1)4. 15(1 and Eg. H4.1|) . 

If we set w{f3) = 5{f3 — (3o), < /3o < 1, we arrive at the "normal" form 
of a fractional diffusion equation, Ea. ()2.5() . with (i = {Sq and Kp = Kt^~I^° . 

Let us prove that the solution of Eg. H4.15() is a PDF. Applying the 
Laplace- and Fourier-transforms in succession, we get: 
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where 

'1 



Ic{st) = / dpisrfwiP). (4.17) 
Jo 

We rewrite Eg. 1)4. 16(1 in the form analogous to Ea. H4.4|) . 

f(k,s) = — due-^'V^^^^^^] = due-''^'^^Gc{u,s) (4.18) 
s Jo Jo 



where 

G^(^,s) = ^^^e-"^^(-) (4.19) 



s 

is the Laplace transform of a function Gc{u, t) whose properties will be spec- 
ified below. Now, p{x,t) can be written in the form analogous to Ea. (|4.6p : 

Pix,t)= r^e-^'^f ^e^' rdue-''^-Gciu,s) = 
7_oo 27r Jsr 27rz Jq 

oo ^—x^/4uKt 



du — Gc{u,t) (4.20) 

v47rnAT 

Similar to the RL case, the function Gc{u,t) is the PDF providing the 
subordination transformation, from time scale t to time scale u. Indeed, at 
first we note that Gc{u,t) is normalized with respect to u for any t. Using 
Eg. lflTT^ we get 

duGc{u,t) = L-^ du^e-''^^^^ =L-i|^| = 1. (4.21) 

To prove the positivity of Gc{u, t) it is enough to show that its Laplace 
transform is completely monotone on the positive real axis ^Sj. The last 
statement is proved by noting that it is a product of two completely mono- 
tonic functions, Ic/s and exp(— n/c). We again demonstrate this for the 
particular choice of w{f3), see Eg. 1)4. 8(1 . For this choice we obtain from 
Eq.dHTTj) 

/c(s) = 6is^i +62S^^ (4.22) 

where bi = Bit^^, 62 = B2T^'-^. Then, Ic/s is completely monotone as a 
sum of the negative powers of s, and Ic itself is a positive function with a 
completely monotone derivative. Thus, exp(—ulc) is also completely mono- 
tone. Thus, we have proved that is completely monotone, and Giii{u,t) is 
a PDF, according to the Bernstein theorem. 

We are interested in the behavior of the MSD, i.e. in second moment of 
the PDF. Using Eg. (|l?T^ . we have 
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Inserting Eq. (|4.22|) into Eq. (|4.2:-{|) one obtains: 

<^^«> - { JIm^tVm^} - f^^- { feT55^} • 

(4.24) 

Recalling the Laplace transform of the generalized Mittag-Leffler function 
Efj,,u{z), > 0, u > 0, which can be conveniently written as jllj 

L {t'^-^E.^^i-Xtn} = Res > \X\'/^ , (4.25) 

we get from Ea. ()4.24() : 

(xHt)) = ^t^^i?,,-,,,,,^-! [-^^/^-'') ■ (4.26) 

To obtain asymptotics at small t, we use an expansion, which is, in fact 
the definition of E^^^{z), see [TS|, Ch.XVIII, Eq.(19): 

°° n 
n=0 ^ ' 

which yields in the main order for the MSD 

J!l"J 'X (4.28) 



For large t we use the following expansion valid on the real negative axis, 
see Uni, Ch.XVIII, Eq.(21): 

^ -n 

which yields 

Since (3i < (32^ we have retarded subdiffusion. It is worthwhile to note that 
the distributed-order equation (|4.15jl with the Caputo derivative always 
describes retarding, or slowing-down, sub-diffusive processes. Indeed, it is 
clearly seen from Eas. (|0^ . that, for w{l3) given by Ea. (|0|) . for 

large s (short times) I{st) behaves as s^^, so that, due to the Tauberian 
theorem, the MSD behaves as t^^, while at long times (small s) I{st) oc s^^ 
and, respectively, (x^) oc t^^ , f3i < ^2- Thus, Ea. (|4.15|) cannot describe the 
accelerating sub-diffusive process. The same arguments applied to Eg. 1)4. 1(1 
with the Riemann - Liouville derivative demonstrate that Eg. 1)4. 1(1 cannot 
describe retarded sub-diffusive processes. 
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5. Distributed-Order Space Fractional Diffusion Equation 
5.1. The "normal" form 

We write the distributed-order space fractional diffusion equation in the 
"normal" form as 

^ = £daK{a)^, p{x,0)=6{x) (5.1) 

where K{a) is a (dimensional) function of the order of the derivative a, 
and the Riesz space fractional derivative d°'/d\x\°' is given by its Fourier 
transform, see Eq. (|:-{.6j) . Setting K{a) = KaQ6{a — ao) we arrive at Eq. llH.ljl . 
In the general case K{a) can be represented as 

K{a) = r-^Kw{a), (5.2) 

where I and K are dimensional positive constants, [/] = cm, [K] = cm^/sec 
and w is a dimensionless non-negative function of a. The equation for the 
characteristic function of Eq. (|5.1|) has the solution 

f{k,t) = expl^-^ daw{a){\k\l)°Y (5-3) 
Note that the normalization condition, 

dxp{x, t) = f{k = 0, t) = 1 (5.4) 



/ 

J — ( 



is fulfilled. 

Let us consider the simple particular case, 

w{a) = Ai5{a — ai) + A25{a — 02), (5-5) 

where < ai < 02 < 2, ^1 > 0, > 0. Inserting Ea. ((^ into Ea. ((^ 
we have 

f{k, t) = exp {-ai \k\'^^ t - aa \k\'^^ t} (5.6) 

where oi = AiK /t^^"-'^ , 02 = A2K /t^^"-"^ . The characteristic function 1)5. 6() 
is a product of two characteristic functions of Levy stable PDFs with the 
Levy indexes ai, Q2, and the scale parameters a]/°'^ and ag^"^, respectively. 
Therefore, the inverse Fourier transformation of Ea. H5.6|l gives the PDF 
which is a convolution of two stable PDFs, 



— ^ — ^ 1 i_ f°° 




02,1 



,(«2i)°2 / 

(5.7) 
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where L^^i is the PDF of a symmetric Levy stable law given by its charac- 
teristic function 

L,,i(A:) =exp(-|A:r). (5.8) 

The PDF given by Ea. ()5.7|) is, obviously, positive, as the convolution of two 
positive PDFs. 

The PDF will be also positive, if the function A[a) is represented as a 
sum of delta-functions multiplied by positive constants, is a positive 
integer. Moreover, if A[a) is a continuous positive function, then discretiz- 
ing the integral in Ea. ()5.1|) by a Riemann sum and passing to the limit we 
can also conclude on the positivity of the PDF in the general case. 

Let us consider the q-ih. moment, Eq. (|^-i.l()|) . at small and large t. After 
inserting Eq. (|5.6j) into Eq. (|H.l()|) we get 

dxl") = ^r(l + q) sin dk (l - e-**W) k^'^-^ , (5.9) 

where 

(t){k) = aik''^ + a2k''\ (5.10) 

In order to get the q-ih. moment at small t, we expand exp (— ai t) 
in power series with subsequent integration over k. As the result we have 
the following expansion valid at g < qi: 

dxl'^) = — (ast)"/"^ sin (^) r(l + q)T f 1 - — x (5.11) 

\ r(i-^)ir; V "2 y j 

for t — > 0. We recall that the radius R of convergence of a series X^^q '^nx'^ 
is determined from the equation 

1/R = lim Icnl"*"^" . 

n— >oo 

Then, by using the Stirling formula, one can see that in case of Ea. H5.11() 
|cn|^''" ~ ^ for n oo, and therefore expansion ()5.1H) is valid 

at all t. The leading term of a series (|5.11() gives for the characteristic 
displacement at small t, 

(^x ~ oc ti/"^ (5.12) 

In order to get the q-th moment at large t we integrate by parts the right 
hand side of Eq. 1)5. 9|) : 

(|xn = ^r(l + g)sin(^)^J(t), (5.13) 
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where 

f'OO 

J{t) = / dkQ{k)e-^^^''^ (5.14) 

JO 

with 

Since at smaU k one has </>(A;) ~ oi/c"^, C(^) ~ oiaiA;"^^'^^"'^, then at large t 
we have exactly the case of Laplace asymptotic integral, see ^7j) Chapter 
III, Theorem 7.1. For J{t) one immediately obtains J{t) = af "'i'?/°i"^r(l- 
q/ai) and 

(Ix^) ~ — (ait)"/-^ sin (^) r(l + g)r f 1 - g < ai. (5.15) 

For the characteristic displacement at large t we get 

~ oc t^/^^ (5.16) 

Therefore, at small times the characteristic displacement grows as t^/"^, 
whereas at large times it grows as Since ai < one encounters an 

accelerated superdiffusion. 



5.2. The "modified" form 

Let us now consider the conjugated form of the equation, namely one 
with an additional distributed-order fractional operator in the l.h.s.: 



daw{a)l " 







d \x 



2-a 



dt 



-K^, p{x,0) = 6{x), (5.17) 



where /, K and w have the same meaning as in Ea. H5.2p . Setting w{a) = 
6{a — ao) we arrive at Ea. (|3.8|) where Ka = = Kl°°~'^. The equation 
for the characteristic function of the solution of Ea. H5.17p reads: 



fik,t) = exp 



Kt/f 



Jg daw{a){\k\ ly 



(5.18) 



(compare with Ea. (|5.3|) '). Note that the normalization condition, Ea. H5.4|) . 
is fulfilled. 

As we did it throughout the present article, let us consider a particular 
case of 

w{a) = Ai6{a — ai) + A26{a — 02), (5.19) 
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< «! < Q!2 < 2. From Eg . (15.181) we get for the characteristic function, 

t 



f{k,t) = exp 



ai 



ifcr 



- + 



a2 



(5.20) 



where ai^2 = Ai^2l'^~°'^-^ /K. The proof of non-negativity of the PDF given 
by inverse Fourier transform of Eq. (|5.2U|) follows along the same lines as for 
the accelerating subdiffusion case. 

The q-th moment is given by Ea. H5.9|) . where 



(5.21) 



We insert Eq. (|5.2H) into Eq. (|5.9p and pass to a new variable ^ = tk°'^ /ai. 
For the integral over k we get 



ax) ai 



(5.22) 



dee 



1 — exp 



a2 0-2 -I 



1 + £2£^-^ 
ai 



At small t we can neglect the term with t in the denominator in square 
brackets. Therefore, 



{\xn=Ciq,ai) 



q < ai, 



where C{q,ai) is given by Ea. (|3.13() . and 



(5.23) 



(5.24) 



for t 0. 

In order to get the q-th moment at large t, we again use the Laplace 
method. Turn to Ea. (|5.13() . where <p{k) is given by Ea. H5.21|) and 



Cik) 



(t)'{k) _ 1 aiaifc-°i +a2a2fc~°^ 
~k^ ~ k^+1 {aik-^^ +a2k-'^2f 



(5.25) 



Since at small k we have (/'(fc) ~ /c"^, C{k) ~ (0:2/02) then at large 

t we again have exactly the case of Laplace asymptotic integral jl7j . For 

J{t) we have immediately J{t) = a2'^^'^^t'^^'^^~^T{l — qja^)-, and 



(|xr)«C(g,a2)(-j 



q < ai, 



(5.26) 
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where C{q,a) is given by Eq. (|;-i.l,Sj) . and 

5x ~ dxl'?)^/'? oc (5_27) 

for t — > oo. 

Therefore, at small times the characteristic displacement grows as t^/"^ , 
whereas at large times it grows as t^^"^. Since ai < 02, we encounter 
retarding superdiffusion. 

6. Conclusions 

Distributed-order diffusion equations generalize the approach based on 
fractional diffusion equations to processes lacking temporal scaling. The 
typical forms of fractional diffusion equations can be classified with re- 
spect to the position of the fractional derivative in time/coordinate in- 
stead of or in addition to the first and second derivatives in the classical 
Fick's form. In the present paper we considered the corresponding forms of 
distributed-order fractional diffusion equations and elucidate the effects de- 
scribed by different classes of such equations. We show that equations with 
the distributed-order fractional operator replacing the corresponding whole- 
number derivative describe processes getting more anomalous in course of 
the time, i.e. the accelerating superdiffusion or retarded subdiffusion. On 
the opposite, equations with additional fractional operators on the "wrong" 
side of the diffusion equation describe processes getting less anomalous, i.e. 
retarded superdiffusion and accelerating subdiffusion. 
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